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We consider an all m-fiber optical modulator based on a ring resonator configuration. The case 
of adiabatic to nonadiabatic transition is considered, wliere tiie geometrical (Berry) phase acquired 
in a round trip along the ring changes abruptly by tt. Degradation of the responsivity of the 
modulator due to finite linewidth of the optical input is discussed. We show that the responsivity 
of the proposed modulator can be significantly enhanced with optimum design and compare with 
other configurations. 
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I. INTRODUCTION 

Optical modulators are devices of great importance for 
optical communication and other fields. In these de- 
vices some external perturbation, e.g. electric or mag- 
netic fields, is employed to modulate the transmission T 
between the input and output optical ports. One of the 
key property of an optical modulator is the responsivity, 
namely the dependence of T on the applied perturba- 
tion. Enhancing the responsivity is highly desirable in 
many applications. As is shown in Ref. 0, the linearity 
of optical modulators imposes in general an upper bound 
on their responsivity. One way of achieving high respon- 
sivity is by employing a resonator configuration with high 
quality factor Q. The multiple back and forth reflections 
occurring in a resonator allow enhancing the responsivity 
in comparison with the case of reflectionless optical path. 
Such a ring resonator was considered recently by Yariv 
0, HI and implemented experimentally Q, Q- It was 
shown that high enhancement is achieved when critical 
coupling occurs, namely when the power entering the res- 
onator from the input port equals the output dissipation 
power. On the other hand, one of the drawbacks of a 
resonator configuration is the limited optical bandwidth. 
In some cases finite linewidth of the optical input Aw 
may lead to broadening of the resonance and thus reduc- 
ing the responsivity. Such broadening can be avoided 
only when Alu/lo << X/QL, where A is the wavelength 
and L is a characteristic length of the resonator. 

In this paper we consider a ring resonator similar to 
the one discussed in 0, 01 ■ However while Ref. 0, 
(3 considered the case of polarization independent evo- 
lution, here we study the case of finite birefringence k (s) 
along the optical path (s is a coordinate along the optical 
path). We first consider the case of adiabatic evolution, 
when K changes slowly. In this case it is convenient to 
express the state of polarization (SOP) in the basis of 
local eigenvectors. In this basis the equations of motion 
of both polarization amplitudes can be decoupled to the 
lowest order in the adiabatic expansion. Next we con- 
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sider the case of adiabatic breakdown, namely the tran- 
sition into the regime where the adiabatic approximation 
does not hold. In this case the geometrical (Berry) 
phase acquired in a round trip along the ring changes 
abruptly by tt. We show that this abrupt change can 
be employed for achieving high responsivity. Note that 
similar adiabatic breakdown was considered in Ref. 
for the case of spin 1/2 electrons in coherent mesoscopic 
conductors with spin-orbit interaction (see also Ref. j^])- 
Such an optical modulator based on adiabatic break- 
down can be implemented in a variety of different config- 
urations. Here we demonstrate these effects by consid- 
ering a relatively simple example of a modulator based 
on a fiber ring resonator having both intrinsic and exter- 
nally applied birefringence. The intrinsic birefringence 
along the ring in our example is linear. As we discuss 
below, it can be induced using a standard polarization 
maintaining fiber being twisting and tapered to realized 
the desired birefringence. The externally applied bire- 
fringence used for modulation is based in our example on 
magneto-optic effect 9], This effect allows induc- 

ing circular birefringence in the fiber, being proportional 
to the Verdet constant characterizing the material and to 
the component of the applied magnetic field along the di- 
rection of propagation. We employ both analytical and 
numerical calculations to study the responsivity of the 
system. We find enhanced responsivity when operating 
in the adiabatic breakdown regime. 



II. FIBER RING RESONATOR 

Consider a fiber ring resonator as seen in Fig. ^ It 
consists of a fiber ring coupled to input and output ports 
using a directional coupler. 

The SOP at each point along the fiber is described as 
a spinor with two components associated with the am- 
plitudes of two orthonormal polarization states. As is 
discussed in appendix A and appendix B, we use the 
local eigenvectors as a basis to express the SOP. The 
associated amplitudes are i?| and E;^ respectively. The 
directional coupler is assumed to have coupling constants 
independent of the SOP. Moreover, the coupling is as- 
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Note that if M is unitary (namely, M ^ = M'') and 
13 holds then^ as expected, S is unitary as well. Note 
also that if M is diagonal (namely, M12 = M21 = 0) the 
following holds 
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t-M22 
1-M22f 



(7) 



FIG. 1: The fiber ring resonator. 



sumed lossless, thus the coupling matrix is unitary 



where 




T 



(1) 



\t\' + \r\'^l. 



(2) 



Integrating the equation of motion along the ring leads 
in general to a linear relation between the amplitudes at 
both ends 
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where 



M 



Mil M12 
M21 M22 



(4) 



Using ^ 13 and |2] one can find a linear relation be- 
tween the amplitudes in the input and output ports of 
the modulator 



£■61 
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T 

E<^i 



(5) 



where the matrix S is given by 



To find the matrix M one has to integrate the equation 
of motion lXol along the close curve defined by the ring. In 
the adiabatic limit, to be discussed in the next section, 
the solution can be found analytically. In the follow- 
ing section the case of adiabatic breakdown is discussed, 
where both analytical approximations and numerical cal- 
culations are employed to integrate the equation of mo- 
tion EHl 



III. THE ADIABATIC CASE 

In the case where the adiabatic approximation can be 
applied the matrix M is given by 



M 



exp(z5|) 
exp{idi) 



(8) 



where 5^ and are given by equations 113231 and \B24\ 
respectively. 

In the more general case the ring may have internal 
loss. Assuming the loss is polarization independent, one 
has 



cxp(i(5|) \ , . 

exp(z,5^)j' 



where < < 1 is real. Thus using [7| 



Ebi 



t- (l-$;)exp(i(5^) 



(10) 



where a E {Tii}- Using the notation t = 
(1 — ^^) exp (i^t), where < < 1 is real, and = 
Sa — Ot one gets 



El' . l-ec-(l-ez)exp(»^) 

^ = -ejexp(..)- (^1) 

Near resonance t9 << 1. Moreover, assuming << 1 
and << 1, one finds 



1 - tM-^ 



t* 



(6) 



|^.exp(.,)^-^-- 



E'' 



id 



(12) 
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Critical coupling occurs when = = ^. In this case 
the transmission amplitude E^^ /E^^ vanishes at reso- 
nance. The transmission probability in this case is given 
by 



1 + [Qdy 



(13) 



where Q — 1/2^. Thus, high responsivity can be 
achieved when operating close to a resonance with high 
Q factor. 



IV. BROADENING DUE TO FINITE 
LINEWIDTH 

As was discussed in the previous section, relatively 
high responsivity can be achieved when operating close 
to a resonance. However, as we discuss below, the price 
one has to pay for that is limited bandwidth. 

Consider the case where the optical input has some 
finite linewidth Acj. As a result the phase factor d will 
acquire a linewidth given by 



M ^ 27r -. 

uj A 



(14) 



Consider the case of a polychromatic optical input and 
assume that the probability distribution of •& is Lorenzian 
with a characteristic width At5 



(15) 



Averaging using this distribution and Eq. 1131 and em- 
ploying the residue theorem for evaluating the integral 
one finds 



T(i9) = / dd'f{d')T{d') 



(16) 



= 1- 



1 



1 



( V 

\l+QAi} ) 



Thus, for this case broadening can be avoided only if 
QAd « 1 or Alj/u; « X/QL. 



V. ADIABATIC BREAKDOWN 

While in the previous case both adiabatic SOP are 
effectively decoupled, we consider now the transition be- 
tween adiabatic and non-adiabatic regimes. 

The birefringence along the fiber ring is described by 
the vector k (s) (see appendix A) . Consider the case 
where in some section of the ring k,{s) is close to the 



degeneracy point at the origin k. = 0. In this case small 
perturbation applied to k (s) can result in a large change 
in the geometrical phase IB20I and IB21I This can be 
seen by considering, for example, the case of a planar 
curve k(s). In this case the solid angle is given by 
— 27rr7,, where n is the winding number of the curve 
K (s) around the origin. As the curve n (s) crosses the 
origin at some point, n changes abruptly by one, leading 
thus to an abrupt change in the geometrical phase. Note 
however that near this transition when |/«(s)| is small 
the adiabatic approximation breaks down and alternative 
approaches are needed. 

As an example for such a transition we consider a ring 
resonator for which the close curve k (s) has the shape 
seen in Fig. 111(c) in the unperturbed case. This curve is 
made of 'half circle' section in the 1-3 plane (the linear 
birefringence plane) and a 'diameter' section along the 
K3 axis crossing the origin. Such a structure can be 
realized by using a polarization maintaining fiber and 
by employing fiber tapering techniques. The half circle 
section can be made out of a Mobius like ring made of 
the polarization maintaining fiber. After welding the two 
ends of the twisted fiber to form the Mobius structure one 
can employ tapering techniques to form the 'diameter' 
section. 

The curve k (s) is perturbed by applying a magnetic 
field on part of the 'diameter' section of the fiber ring. 
Such a perturbation contributes circular birefringence in 
the K2 direction (see Fig. 13 (a) and (e). The rela- 
tively high value of the Verdet constant in common op- 
tical fibers allows significant magneto-optic effect with 
moderate applied magnetic fields. While the adiabatic 
approximation totally breaks down in the unperturbed 
case of Fig. |5|(c) when the curve k (s) crosses the origin, 
the perturbation transforms the system into the regime 
where adiabaticity holds. As is shown below, the respon- 
sivity of the system is relatively high when operating near 
this transition between the adiabatic and non-adiabatic 
regimes. 

The 'half circle' section is analyzed in appendix C. As 
can be seen in Fig. |S1 the Zener transition probability 
Pz vanishes for a series of points denoted as A„. In our 
example we chose A to be the first zero of Pz (A), namely 
A = Ai = 1.022. One advantage of choosing one of the 
zeros oi Pz (A), where Pz obtains a local minimum, is the 
fact that Pz is only weakly affected by small deviations of 
K (s) from the ideal 'half circle' curve. For the parameter 
7 we chose the value 7=1. As can be seen from Fig. 
121 (d) for this choice the evolution along the 'half circle' 
section transform the polarization vector on the Bloch 
sphere from the pole on the negative Pz axis to the op- 
posite pole on the positive Pz axis. The fiber length of 
this section is 2A1/7. 

The rest of the fiber ring has a birefringence given by 
K (s) = Kq (s) -|- Ki (s), where Kq (s) is the unperturbed 
birefringence forming the 'diameter' section and Ki (s) 
is the perturbation induced by the magnetic field. The 
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where A = 50 and B = 0.6 in our numerical example. 
Thus, the magnetic field is applied to a fiber section of 
length 2B(3/"f and drops down to zero abruptly outside 
this section (due to the large value chosen for the pa- 
rameter A). The coupling constants in the numerical 
example are — 10^^ and = 10^''. 

The equation of motion along the fiber ring is inte- 
grated numerically as described in appendix A. This 
allows calculating the evolution of the polarization vec- 
tor on the Bloch sphere (see Fig. [3 (b), (d), and (f)). 
The same calculation yields also the matrix M. The 
off-diagonal matrix elements allow calculating the Zener 
transition probability |Mi2|^ = |M2i|^ (see Fig. 121(a) 
solid line). The curve shows the gradual transition be- 
tween the non-adiabatic limit where \a\ << 1 and the 
adiabatic limit |q;| >> 1. An approximated analyti- 
cal expression for the Zener probability in a similar case 
where the curve n (s) is an infinite straight line was de- 
rived in appendix C. The result in Eq. IC25l can be used 
to estimate approximately the Zener transition probabil- 
ity for the present example 



FIG. 2: The birefringence k (s) and polarization P (s) along 
the fiber ring. Plots (c) and (d) shows the unperturbed case, 
while in (a) and (b) the perturbation parameter is a = —0.2, 
and in (e) and (f) a = 0.2. 



unperturbed part is assumed to be given by 



Pz 



exp 



'Al7 



(19) 



The estimate in Eq. E| is shown in Fig. |31 (a) as a 
dashed line. The deviation between the numerical and 
analytical results is originated mainly by the fact that the 
straight line section in k (s) is finite while the analytical 
analysis assumes an infinite straight line. Moreover, the 
analytical result is expected to hold only in the limit 
where |pz| << 1 as it is evaluated only to lowest order in 
the adiabatic expansion. 

Figure 13 (b) shows the phase of both diagonal matrix 
elements of M. In both cases the phase changes abruptly 
by TT near a = 0. This is originated by the sharp change 
of the solid angle by 27r near a = (see Eq. IB23l and 
IB24|) . The optical modulator discussed in the present 
work employs this sharp change to achieve high respon- 
sivity. 

Figure |3| (c) shows the transmission probability into 
both SOP, Pii = IS-iil^ (solid line) and P21 = \S2i\^ 
(dashed line) of the entire modulator. For both cases, the 
fuU width half maximum (FWHM) is Aa = 5.1 x lO'^. 



Kq (s) 



A172 



(17) 



where \s\ < P/j. In our numerical example the di- 
mensionless parameter /3 is given the value P = 5. The 
perturbation due to the applied magnetic field gives rise 
to birefringence given by 



Ki (s) 



0,- 



1 -I- exp A 



f] -B^ 



,0 



(18) 



VI. DISCUSSION 

As we have seen, the ring resonator can serve as an 
optical modulator with high responsivity when operated 
near one of its resonances. Two regimes of operation 
were considered, the adiabatic one, and the non-adiabatic 
one. In what follows we compare between both regimes 
by considering the following points. 

Optical Source Linewidth - In the adiabatic limit, 
when the equations of motion in the adiabatic basis be- 
come decoupled, the only effect of the external pertur- 
bation is on the phases acquired along the fiber ring. 
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VII. SUMMARY 



In the present work we study an optical modulator 
based on a fiber ring resonator. Both adiabatic and 
non-adiabatic regimes of operation are considered. We 
find that operating close to the point where the geomet- 
rical phase changes abruptly by tt can allow relatively 
high responsivity, even when coupling is not set to be 
critical. Our example deals with a particular configu- 
ration of optical fiber ring with both intrinsic and exter- 
nally applied birefringence. However, the same ideas can 
be implemented with other optical waveguides and other 
birefringence mechanisms. 



FIG. 3: Dependence on the perturbation amplitude a. (a) 
Zener transition probability, calculated numerically (solid 
line) and estimated using Eq. IC25I (b) Phase of M\\ (solid 
line) and of M22 (dashed line), (c) transmission probabil- 
ity into both SOP, Pn = ISn]^ (solid line) and P21 = l&ij^ 
(dashed line). 
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The dependence of the dynamical phase on wavelength 
gives rise to broadening of resonances when operating 
with an optical input having finite linewidth. In the 
general nonadiabatic regime, however, the external per- 
turbation can affect not only the phase factors but also 
the SOP as it evolves along the close fiber ring. The 
later, being wavelength independent gives rise to a mod- 
ified dependence on the optical source linewidth. 

Critical Coupling - In the adiabatic regime full mod- 
ulation between zero and one of the transmission prob- 
ability T is possible only when critical coupling occurs, 
namely t,^ = £,1, (s'^e Eq. 112(1 . In practice, fulfilling this 
condition when = << 1 is difficult. However, this 
condition is not essential in the general non-adiabatic 
case. As can be seen in Fig. O (c) full modulation is 
achieved, even thought for this example = 100^;. 

Responsivity - The responsivity of the ring resonator 
device can be characterized by the FWHM and height 
of the resonance near which the device is being oper- 
ated. As was discussed above, Aa = 5.1 x 10"'^ for 
the example presented in Fig. |3| For the same pa- 
rameters the FWHM of the resonances in the adiabatic 
regime \a\ >> 1 can be evaluated using Eq. 1121 vielding 
Aa — 1.7 X 10^'^. However, as was discussed above, 
since the coupling is not critical, the modulation is not 
full in the adiabatic case. Note that in general the re- 
sponsivity has an upper bound imposed by the linearity 
of the system Qj. It can be shown that for both cases, 
the obtained responsivity is of the same order as the up- 
per bound. A future publication will discuss this point 
in more details. 



APPENDIX A: SOP EVOLUTION ALONG A 
FIBER 

Consider an optical fiber winded in some spacial curve 
in space. Let r (s) be an arc-length parametrization of 
this curve, namely the tangent s =dr/ds is a unit vector. 
The normal unit vector £> and the curvature n are defined 
as ds/ds ~ KU. One can easily show that i> • s = by 
taking the derivative of s • s = 1 with respect to s. The 
vectors s, i> and the binormal unit vector, defined as 
b = s X i>, form a local triplet orthonormal coordinate 
frame known as Serret - Frenet frame [13, [11 (see Fig. 
0J|. By taking the derivative of s • i> = with respect 
to s one finds s-du/ds = —k. Similarly, by taking the 
derivative of b • i> = with respect to s one finds b • 
dCz/ds — —i>- dh/ds. Using the definition b = s x /> one 
finds db/ds = s X (ij>/(is. Thus s-db/ds = 0. Moreover, 
by taking the derivative of b • b = 1 with respect to s 
on finds h-db/ds = 0. Thus dh/ds is parallel to C>. 
The torsion r is defined as dh/ds = —tP. The above 
definitions and relations can be summarized as follows 

The equation of motion along the optical ray defined 
by the fiber can be obtained using the transport equation 
of geometric optics 1^ for the electric field phasor Eg 

2 (VV' • V) Eo + Eo [V V - V (In ^) • VV-] (A2) 
+2 [Eo- V(lnn)]VV' = 0, 



6 



Equation IA7I is known as Rytov's law [3| ■ In the 
more general case where other birefringence mechanisms 
are present the equation of motion reads 

^ |e) = lie |e) , (A9) 

where K, = tCg+lCf, and /C/ is the birefringence in the 
fiber due to intrinsic structure or due to elasto-optic or 
electro-optic of magneto-optic effects. 

In a lossless fiber the matrix JC is Hermitian. For this 
case it is convenient to expresses K, as 



K, — kni 



(AIO) 



where / is the 2 by 2 identity matrix, feg is a real scalar, 
K = K2, K3) is a three-dimensional real vector and the 
components of the Pauli matrix vector cr are given by: 



FIG. 4: The Serret-Frenet frame. 



where ip is the eikonal, 
and /i is the permeability. 
Eo 



n is the index of refraction, 
Define the unit vector eg = 



CTi = 

The s 
motion 



Eo • Eq in the direction of Eq. In terms of §0 the state at some s > so 



'(All) 

evolution operator u (s, so) of the equation of 
9] relates an initial state |e (so)) with a final 



transport equation reads 



ds 



eo = -K(eo • />)§. 



(A3) 



Expressing the unit vector §0 in the Serret - Frenet 
frame 



eo = e^v + e^b, 
one finds using IXT1 



(A4) 



|e (s)) = u (s. So) |e (so)) 
It can be expressed as 



u (s, So) 



lim 



N 

n exp 



(A12) 



(A13) 



where As = s — so, and s„ = sq + nAs/N. For a finite 
N the above expression can be used as a numerical ap- 
proximation of u (s. So). For calculating the exponential 
terms in lAIBl it is useful to employ the following identity 



ds 



s ds V / 



(A5) 



thus, using the Dirac ket notation 



(A6) 



exp (ixIC) = exp {ik^x) [/cos [ax) + ik ■ cr sin (ax)] , 

(A14) 

where the notation of Eq. lAlOl is being used, and k — 
ka where k is a unit vector and a = \k\. 

The normalized SOP |e) can be represented as a point 
on the Bloch sphere indicating the expectation value of 
the Pauli spin vector matrix, namely 



one finds 



= e (T e 



(AI5) 



— |e) ilCg |e) , 



(A7) 



where the geometrical birefringence ICg is given by 



APPENDIX B: THE ADIABATIC CASE 

To establish notation we review below the main results 
of Rcf . G] . Consider the differential equation 



(A8) 



ds 



(Bl) 
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where \'ip) represents N dimensional column vector and 
K- = K. (s) is N X N Hermitian matrix. For any given 
value of s the Hermitian matrix IC (s) has a set of or- 
thonormal eigenvectors 

IC\nis))=Kn{s)\n{.s)), (B2) 
where n — 1,2, N and 

(n (s) |m (s)) = Snm- (B3) 
The solution can be expanded as follows 



L Jo 



(B4) 



Substituting in lBll vields 



n 

(B5) 

where upper-dot represents derivative with respect to 
s. The off-diagonal terms, given by 



(m (s) |ri (s)) 



(m(g)|/C|n(s)) 

i^n (s) - Km 



(B6) 



where m ^ n, are neglected in the adiabatic approx- 
imation. The resulting decoupled set of equations are 
easily solved 



(B7) 



(B8) 



Consider now the two dimensional case N ^ 2. Using 
the notation of Eq. lAlOl and the notation k — ka, where 
/« is a unit vector, given in spherical coordinates by 



a„ (s) = a„i (0) exp (z7„) . 
where the real phase 7^ is given by 

s 

7m =M ds' {m (s') \m {s')) ^ 



and the following holds (T I T) = (i U) = 1, (T U) 
0, and 



/C|T) = (fco + a)|T) 



/C li) = (fco - a) li) 



(B12) 



(B13) 



The eigenstates \n) (where n E {?> i}) are independent 
of ko, thus: 

7,„ = i J ds (m (s) \rji (s)) = i J dK ■ (m \m (k)) . 

Sl Ki 

(B14) 

Using the expression for a gradient in spherical coor- 
dinates one finds 



(T|v.|T> = -^ctg( 



aiv.ii) = ^ctg0. 



(B15) 



(B16) 



For the case of a close path. Stock's theorem can be 
used to express the integral in terms of a surface integral 
over the surface bounded by the close curve k (s) 



7„j ^ i <j) dn- (m| Vk \m) ~ i j da • (V x (m| \'m)) . 

(B17) 

Expressing the curl operator in spherical coordinates 
one finds 



Vx(T|v.|T) = ^^, 



vxaiv,u) = - 



I K 



(B18) 



(B19) 



a — (cos if sin 9, sin ip sin 0, cos ( 



finds 



IC = fco/ + a 



cos9 sin exp (— zi/?) 



sin0exp(i(p) — cos0 
The orthonormal eigenvectors are chosen to be 



(B9) 



(BIO) 



exp(-f) 

■ ) 
(Bll) 



cos J exp ( ^ 



and 



7t 



1 /■ K 1 

- da - — 5- = — 
2; 2 



da • — 5- = 

2 



(B20) 



(B21) 



where is the solid angle subtended by the close path 
K (s) as seen from the origin. Due to the geometrical 
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nature of the last result, the phase factors 7^ and 7^ are 
called geometrical phases. 
Thus 



where 



5^ = --+J^ ds'[ko{s')+a{s')] 



(B23) 



Si = ^+ ds'[ko{s')-a{s')], (B24) 



APPENDIX C: ZENER TRANSITIONS 



The set of equations IB5I for the two dimensional case 
N = 2 can be written in a matrix form as follows 



ds V "-l 



(CI) 



T 



-cxp(-i/3) ( i 



-exp(z/3)(t li 



by 



Thus, to lowest order the transition probability is given 



Pz 



ds'exp[z (/3 + 7,)] (t li) 



(C6) 



Consider the case where k = 
a (cos sin f?, sin sin 0, cos 0) is planar with ip — ccmst. 

Using En] 



thus 



Similarly 



(tli) 



= 0, 



thus 7^ = 0. Using the above results 



d9 exp {iQ 



(C7) 



(C8) 



(C9) 



(CIO) 



where 



where 



/o 

= -2 



ds' [K,{s')-K, {s')\ 
ds' \n\ . 



(C2) 



In the adiabatic limit the off-diagonal matrix elements 
are considered negligibly small, and consequently no 
transitions between the adiabatic states occur. To calcu- 
late the transition probability to lowest order we consider 
the off diagonal elements as a perturbation 0. The so- 
lution of the unperturbed problem is given by 



Oi (s) = (O)exp {i-i^ 



AT (0) exp [ill) 



(C3) 



(C4) 



Assuming that at some initial point sq the system was 
in the ||) state, we wish to calculate the probability to 
find the system in the |t) state at s > sq. Lowest order 
correction is obtained by substituting the unperturbed 
solution in lCll 



a| = —a^ (0) exp 



/3 + 7, 



(C5) 



m 

-2 I ds'\K\ 
10 



(CII) 



1. The case where k (s) is half circle 

Consider the case where /C = k • cr, where 

'^(s) = 7(\/a2-(7s)',0,7s) , (C12) 

7 is a non-negative real constant with dimensionality 
of 1/length, A is a non-negative dimensionless real pa- 
rameter, and |s| < A/7. 

The Zener transition probability is calculated for the 
case A > 1 to lowest order in the adiabatic expansion. 
The following holds 



7s 

X' 



(C13) 



and \k\ — 7A, thus 



MS) 

C (9) = -2 ds' \k\ = -2A2 cos 0, 
Jo 



(C14) 
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FIG. 5: Example of numerical integration of the equation 
of motion for the case A = 5. On the left the curve k; (s) 
is shown and on the right the evolution of the polarization 
vector p (s) on the Bloch sphere is seen. 



and 



/ d0 exp (— 2iA^ cos ( 



Using the identity 



/ d9 exp {iz cos 6) — nJo {z) , 





one finds 



(for A > 1) 



(C15) 



(C16) 



(C17) 



Fig. |S1 shows an example of numerical integration of 
the equation of motion for the case A = 5. Fig. |S1 (a) 
shows the 'half circle' k (s) curve and Fig. [51(b) shows the 
evolution of the polarization vector on the Bloch sphere. 
FigureElshows a numerical calculation of the Zener tran- 
sition probability as a function of the parameter A. As 
can be seen in Fig. |B| Pz vanishes for a series of points 
we denote as A„ (n — 1, 2, 3...). The first zero of p^ is at 
Ai — 1.022. Note however that even though p^ vanishes 
at the points A„, the evolution becomes truly adiabatic 
only when A >> 1. 

Comparing Eq. IC17I with the numerical solution seen 
in Fig. shows, as expected, good agreement for A > 1. 
For the range < A < 1, however, we find that the 
following can serve as a good approximation 



Jr 



/7rA2 



(for < A < 1) . 



(C18) 



FIG. 6: Zener probability vs. the parameter A calculated 
numerically. 



2. The case where k (s) is a straight line 

We calculate below pz for the case K, — k ■ tr, where 
K (s) is a straight line 



k(s) = A(0,l,7s) 



(C19) 



where A and 7 are real constants independent of s. 
For the present case one has 



C(5) 



2A f ds'Jl + {-fs'f (C20) 

JO 

^s\J 1 + (75)^ + sinh^"'" (7s) 



and 



dO exp {iC,) 
1 



cosh 2; 



exp 



-I — — smh 2z + z 
7 V2 



(C21) 



dz. 



In the limit A/7 —s- 00 the phase oscillates rapidly and 
consequently —t 0. The stationary phase points z„ in 
the complex plane are found from the condition 



= — I isinh2z + z ) =cosh2z + l, (C22) 
dz \ 2 ' 



thus 



(C23) 
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where n integer. Note, however that the term 1 / cosh z 
has poles at the same points. Using the Cauchy's theo- 
rem the path of integration can bo deformed to pass close 
to the point Z-i = —1^/2. Since the pole at Z-i is a 
simple one, the principle value of the integral exists. To 
avoid passing through the pole at z-i a trajectory form- 
ing a half circle " above" the pole with radius e is chosen 
were e ^ 0. This section gives the dominant contribu- 
tion which is inR, where R is the residue at the pole. 
Thus one finds: 



The prcfactor in front of the exponent is determined 
by requiring = 1 in the limit A << 7, thus 



Pz exp I —TT 



A 



(C25) 



2 7. 



d9 exp 



~ exp - 



TT A 



(C24) 
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